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The X-ray luminosity function of AGN at 2: ^ 3 
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ABSTRACT 

We combine Lyman-break colour selection with ultradeep (> 200 ks) Chandra X-ray 
imaging over a survey area of ~ 0.35 deg^ to select high redshift AGN. Applying 
careful corrections for both the optical and X-ray selection functions, the data allow 
us to make the most accurate determination to date of the faint end of the X-ray 
luminosity function (XLF) at z ^ 3. Our methodology recovers a number density of X- 
ray sources at this redshift which is at least as high as previous surveys, demonstrating 
that it is an effective way of selecting high z AGN. Comparing to results at z = 1, 
we find no evidence that the faint slope of the XLF flattens at high z, but we do find 
significant (factor ^ 3.6) negative evolution of the space density of low luminosity 
AGN. Combining with bright end data from very wide surveys we also see marginal 
evidence for continued positive evolution of the characteristic break luminosity L». 
Our data therefore support models of luminosity-dependent density evolution between 
z = 1 and z = 3. A sharp upturn in the the XLF is seen at the very lowest luminosities 
(Lx < lO^^'^erg s~^), most likely due to the contribution of pure X-ray starburst 
galaxies at very faint fluxes. 
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1 INTRODUCTION 

' Understanding the evolution of active galactic nuclei (AGN) , 
and thus the history of accretion in the universe, is an im- 
■ portant problem in astrophysics. Accurate measurements to 
' determine the shape of the luminosity function of AGN, and 
how this evolves with redshift, provide the key data to probe 
the changing distribution of AGN activity, and thus the im- 
portance of accretion onto super-massive black-holes at vari- 
ous times throughout the history of the universe. Such inves- 
tigations require large, unbiased samples of AGN spanning 
a range of redshifts and luminosities. 

X-ray surveys are highly efficient at selecting AGN 
over a large range of redshifts (e .g. iBarger et al. l2003l : 
lEckart et al]|2006l : iBrusa et al.ll2007l ). inc luding uno b scure d 
and moderately obscured AGN (e.g. iGilli et al.1 |2007| ). 
and very faint AGN where the light of the host galaxy 
overwhelms any optic al signature (e.g. iMoran et al. |2002| : 
ISevergnini et al.ll2003h . Thus X -ray selected samples are rel- 
atively unbiased against all but the most heavily obscured 
objects. Large efforts have been made to perform follow- 
up spectroscopic observations of X-ray sources detected in 
various surveys and measure their redshifts, allowing the 
evolution of the X-ray Luminosity Function (XLF) to be 
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investigated. AGN are found to be a strongly evolving pop- 
ulation, with the number density decreasing dramatically 
from z 1 to the present day , similar to the star formation 
rate (e.g. I Merloni et al.|[2004 ) . Thus accretion activity, par- 
ticularly in more luminous systems, was considerably more 
prevalent when the universe was about half its current age. 
This behavio ur can be describe d by Pure Luminosity Evo- 
lution (PLE. lBarger et al.|[2005l '). in which the XLF retains 
the same shape (a double power-law with a break at a char- 
acteristic luminosity, L*), but shifts to higher luminosities 
as redshift increases. The most recent data indicate a com- 
plex picture, in which the redshift at whi ch the AGN space 
density peaks is a function of luminosity l|Ueda et al.ll2003l : 
iBarger et al] l2005l : iHasinger et al.ll2005h . A number of au- 
thors have proposed that a Luminosity Dependent Density 
Evolution (LDDE) p arameterization is n e cessary to describe 
such evolution (e.g . iMivaii et al.l l2000l: lUeda et al.1 l2003l : 
iHasinger et al. l l2005l : lLa Franca et aLlbOOa T In this scheme, 
the shape of the XLF changes with redshift, and various 
theo retical models have attempted to explain such evolutio n 
(e.g. iMerlonilliooi : iHopkins et aLll2006l : iBabic et al.ll2007l ). 



remains unclear. 

10"*^ erg s"^) 



The behaviour at higher redshifts 
While the most luminous AGN (Lx ^ lU'"erg s ' ) ap- 
pear to conti nue to increase in number density at very hig h 
redshift (e.g. IHasinger et al. I l2005l : ISilverman et~aLl |2005| ). 
very little is known about the behaviour of lower luminos- 
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ity AGN. Very deep X-ray data are required to probe lower 
luminosities at these high redshifts, and thus there are few 
suitable fields and the sample sizes are small. In addition, 
the optical counterparts of faint X-ray sources are generally 
extremely faint, making spectroscopic identification difficult 
or sometimes impossible with current instrumentation. Even 
the most intensively-observed fields are therefore highly in- 
complete spectroscopically. In the Chandra Deep Fields - 
North and -South (CDF-N; CDF-S), the deepest Chandra 
observations to date with the most intensive p rogrammes of 
spect roscopic follow-up, 87 pe r cent (CDF-N : Barger et al.l 
I2OO3I ) and 78 per cent (CDF-S: Is^kolv et al.ll2004l ) of X-rav 
sources with i? < 24 have spectroscopic redshifts. This com- 
pleteness falls rapidly for objects with fainter optical coun- 
terparts, however, giving a total of 56 per cent and 39 per 
cent spectroscopically identified X-ray sources in the entire 
CDF-N and CDF-S samples respectively. This can severely 
bias deter minations of the X LF, particularly at high red- 
shifts fe.g. lBarger et al.ll2005l ). The key issue in determining 
the high redshift evolution is therefore how one deals with 
the sources which are too faint to identify. 

Two general approaches have been applied to address 
this. The first is to try to maximise the redshift complete- 
ness of the X-ray detected samples. Generally, this involves 
setting a relatively high X-ray flux limit, where spectro- 
scopic completeness is higher, and augmenti ng the spec- 
troscopy with photometric redshifts (e.g. ,Ueda et al.l 120031 : 
iBarger et aLll2003l : IZheng et al.ll2004) . This allows redshifts 
to be assigned to most of the X-ray objects, and reduces the 
uncertainty in the d etermination of the h i gh-redshift XLF 
llBarger et al.l l2005l : iHasinger et"al] l2005l : iBarger fc Cowid 
I2OO5I ). In the context of our present study, setting a high 



X-ray flux limit is highly undesirable, because the objects 
of interest (i.e. low luminosity AGN at high redshift) are 
at the limits of detectability even in the deepest X-ray sur- 
veys. In addition, there are considerable uncertainties in the 
redshift determinations, particularly for AGN. Finally, it is 
extremely difficult to correct accurately for any residual in- 
completeness in such samples, given that the success and 
failure rates of both spectroscopic and photometric redshifts, 
and errors in the photometric redshifts, depend on many 
complex factors. 

A comp l etely different approach was taken by 



iNandra et al. ([2005a V They used the Lyman-break surveys 
of ISteidel et al. I l|2003D to identify objects in a narrow red- 
shift slice around z ~ 3, based on their broad-band optical 
colours in three filters, with minimal contamination from 
lower redshift sources. These were cross-correlated with X- 
ray detections in 2 deep X-ray fields to identify the pres- 
ence of AGN. As the optical sample is limited to Lyman- 
break selected objects, lower significance detections of X-ray 
sources could be included, without increasing the chance 
of spurious alignments, thus probing to the maximum X- 
ray depth of the available data. A sample of 10 AGN of 
moderate luminosity (Lx = 10'*''~'*^'^erg s~^) were iden- 
tified, 9 of which have been spectroscopically confirmed. 
This colour pre-selection approach samples an incomplete, 
but well defined cosmologi cal volume, which can be cal- 
culat e d from simula ti ons |Steidel et al.l 1 19991 : iHunt et al] 
I2OO4I ). iNandra et aP (|2005af r presented the space density 
of moderate lumi nosity AGN a t z ^ 3, using the survey 
volumes given bv ISteidel etHI |l99i). The result was sie- 



nificantly higher than previous w ork, and prompted the re- 
sponse of lBarger fc Cowi^ l|2005h . using the most up-to-date 
multi-band optical and near infra-red photometric redshifts 
in the CDF-N to reduce incompleteness, increasing their 
result, br inging it to c o nsisten cy with the completeness- 
corrected INandra et al.l l|2005al ) value. In the present pa- 
per we extend and improve on this Lyman-break method 
of determining the high z XLF, using three additional deep 
Chandra fields to increase our numbers of X-ray detected 
Lyman-break galaxies (LBGs), and improving the complete- 
ness corrections. 

In section [2] we describe the optical and X-ray data for 
our 5 fields, including the data reduction, selection of the 
Lyman-break sample and our X-ray source detection proce- 
dure. In section [3] we describe our calculations of the X-ray 
selection function, and the optical selection function, tai- 
lored for our optical data and the selection of AGN counter- 
parts. Section|3]describes our maximum-likelihood fitting to 
determine the XLF, carefully correcting for X-ray and opti- 
cal incompleteness, and our improved method to determine 
binned estimates, accounting for the varying sensitivity of 
our observations. We compare our results to previous work 
in section [5] and discuss the implications for AGN activity 
at high redshift. 

A fiat cosmology with Q.a=Q-7 and h = 0.7 is adopted 
throughout. 



2 DATA/SAMPLE 

The results presented in INandra et all (|2005al l were based 
on the 2 fie l ds fr om the large Lyman-break survey of 
ISteidel et al] (|2003l ) with sufficiently deep X-ray data: the 
Huhhle Deep Field-North (HDF-N), fully contained by the 
Chandra Deep Field-North with ~ 2 Ms exposure time, and 
the Groth-Westphal Strip (GWS) with ~ 200 ks exposure. 
For this work we increase our sample using 3 additional 
fields with ^ 200 ks Chandra imaging with the ACIS-I in- 
stru ment: the Lyii x l|Stern et aLll2002l ) and LALA CETUS 
( Wang et all |2004| ) archival fields, and EGSl (part of the 
AEGIS-X survey, Nandra et al., in preparation), in addition 
to increasing our optical coverage of the HDF-N. 



2.1 Optical data 

2.1.1 Observations and data reduction 

We have obtained deep opt ical imaging of our X-ray fields 
through Un, G and 7^ filters (|Steidel fc HamiltonI 19931 ). suit- 
able for Lyman-break selection at z ~ 3. 

Imaging of the Lynx and LALA CETUS fields was per- 
formed on the William Herschel Telescope (WHT) over 5 
nights in November 2005, using the Prime Focus Imag- 
ing Camera which consists of a mosaic of two 2K x 4K 
EEV CCDs. The 16.2x16.2 arcminute field of view is weU 
matched to the Chandra ACIS-I field of view. The observa- 
tions were performed as individual exposures of 500, 1000 
and 1200 seconds in f/„, G and TZ respectively, with the 
telescope being dithered by ~ 30 arcsec between exposures, 
to compensate for gain variations, bad pixels and the gap 
between the chips. Un band data were always obtained clos- 
est to meridian to minimize the atmospheric extinction. 
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The data were flat-fielded using twilight sky flats. For the 
C/n-band data additional flat-fielding was performed using 
dark sky 'super-fiats', produced by median combining the 
dithered exposures with objects masked out. Fringe removal 
was performed on the 7^-band image using standard tech- 
niques. 

A number of spectrophotometric standard stars from 
the [Oki (|l990l ) catalogue were observed each night and 
were used to calculate photometric zero-points, correct ed for 
Galactic extinction, as described in lSteidel et al.l l|2003l ). The 
images were astrometrically calibrated Jay matching objects 
in the USNO-A2.0 astrometric catalog (|Monet^l998. '). match- 
ing over 100 objects in each field and using a polynomial 
solution to map the focal plane to the catalogue positions. 
The solution was refined for each of the individual exposures 
to allow for small offsets and rotations. The individual ex- 
posures were then resampled to a common pixel scale. Final 
stacks in each band were produced by scaling each frame 
to match the zero-point and magnitudes of objects in the 
photometric exposure observed at the lowest air-mass. This 
allowed nonphotometric data to be included in the stacks, 
although exposures with poor seeing (> 1.5 arcsec) were 
excluded. Finally the stacked images were smoothed with 
a Gaussian filter so the the FWHM of stars matched in 
each photometric band (in practice this meant smoothing 
to match the l7„-band image quality). 

The EGSl field was observed with the Palomar 200 
inch telescope in May 2005, using the Large Format Cam- 
era, which consists of a mosaic of six SITe, back-side illu- 
minated, 2048 X 4096 pixel CCDs. Individual exposures of 
600 (7^), 1200 (G) and 1800 (C/„) seconds were obtained, 
with the telescope dithered inbetween. The data were fiat- 
fielded using dome fiats for the 71 and G band images, and 
dark-sky fiats for the C/„-band. Fringe removal was not nec- 
essary for the 7?,-band with the SITe CCDs. The astrometric 
and photometric calibration, and production of stacks was 
performed as for the WHT data. The EGSl data are of the 
poorest quality in terms of depth and seeing in our sam- 
ple, necessitating the inclusion of poorer seeing data in the 
stacks. 

For the HDF-N field we use the higher qu ality, larger 
area imaging of GOODS-N region, used by iLaird et ahl 
ll2006ll, rath e r tha n the origi nal imaging presente d by 
ISteidel et all l|2003l ) and used in iNandra etall (|2005al ) . De- 
tails of th e observations and d ata reduction have been pre- 
sented bv iReddv et al.l (|2006| ). [/-band optical data were 
obtained with the KPNO Mosaic imager by the GOODS 
team, and transformed to Un magnitudes. G and 7?.-band 
data were obtained on the Keck I telescope with the Low 
Resolution Imaging Spectrograph. 

The optical data for the GWS used in this paper 
are identical to th e 'Westphal' field data presented by 
ISteidel et al] (|2003h . The data were obtained on the Kitt 
Peak 4m Mayall telescope and the Palomar 200 inch in May 
1996 and March 1997 respectively. T he data reduction is 
fully described bv ISteidel erall l|2003l ). and is almost iden- 
tical as our other fields. 

The fields observed, telescopes used, the total expo- 
sure time (after excluding poor seeing data), and the aver- 
age FWHM of stellar objects in the stacked images (before 
smoothing) are summarised in Table [1] 



2.1.2 Source detection, photometry and LBG sample 
selection 

We perform source detection, deblending and photome- 
try f or all of our fields, using the SExtractor software 
(Bcrt in fc Arnouts| [i996V Objects were detected in the TZ- 
band image, after smoothing with a 2 pixel FWHM Gaussian 
kernal, if they contained 5 or more connected pixels with fiux 
exceeding 1.5(trms above the local sky background in the 
smoothed image, where ctrms is the pixel to pixel RMS noise 
in the unsmoothed image. Isophotal apertures were defined 
using the detection criteria. SExtractor also attempts to 
deblend objects with merged isop hotal apertures, but still 
resolvable as separate objects (see lBertin fc ArnoutsI [l996l . 
for details). The same isophotal apertures were then trans- 
ferred to the images in other bands to measure isophotal 
magnitudes, which are used to calculate the C/„-G and G- 
TZ colours. However, the total magnitude we refer to when 
quoting TZ magnitudes is measured within a fiexible aper- 
ture w hich maximises the light collected, as defined bv lKronI 
(|l980h . 

A sample o f z ~ 3 LBGs is se lected using the colour 
criteria given in ISteidel et al] (|2003D : 



G-TZ < 1.2 

Un-G >G-7^-l-1.0 

Additionally objects are required to have 



(1) 
(2) 



19.0 7^ < 



25.5 (Lynx, LALA CETUS) 
25.0 (EGSl) 



and be detected with 3a significance or greater in the G 
band. The bright limit removes sources which may have con- 
tained saturated pixels. The faint limit ensures only high sig- 
nificance detections with accurate photometry are included 
in the LBG sample. We define la limits by A^'pi'^ ctrms where 
Npix is the number of pixels in the isophotal detection aper- 
ture. There is no requirement on the Un band magnitude, 
but objects with a Un band fiux less than the la limit for 
their detection aperture are noted as "un-detected" and as- 
signed a magnitude corresponding to the la limit. The EGSl 
data are of poorer quality and not as deep as the other fields. 
Thus, detection completeness is lower, and the photometry 
is less accurate, which results in more objects being scat- 
tered in or out of the LBG selection box, particularly at 
faint magnitudes. These effects will be accounted for in our 
calculations of the optical selection functions (section 13. 2|) . 
However, we increase the faint limit to TZ< 25.0 in EGSl to 
prevent a high number of faint contaminants being scattered 
into our sample. 

Our detection routine using SExtracto r results in a 
different sample of LBGs to that presented bv lSteidel et al.l 
(l2003f ) for HDF-N and GWS, using FOCAS to perform the 
detection and photometry. The overlap between our new 
samples and those previously presented is ~ 50 per cent in 
the GWS and ~ 40 per cent in the HDF-N(in the area cov- 
ered by both data sets). A number of sources which were pre- 
viously selected as LBGs no longer satisfy the colou r crite- 
ria; m any now satisfy the BX/BM criteria (lAdolbcrg er et al.l 
l2004h . which selects objects in the region of colour space be- 
low the LBG selection, corresponding to objects at slightly 
lower redshifts. A number of additional objects are found 
in our LBG sample (in the same area covered by the orig- 
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inal observations). Such scatter between different data sets 
and photometry is expected, particularly for fainter objects, 
and is accounted for in our calculation of optical selection 
functions (see section [S}. While combining the previous and 
updated samples would increase our number of LBGs, the se- 
lection function would no longer be well-defined. The HDF- 
N and GWS LBG samples used for this work therefore con- 
sist of objects selected by our routines only. 

In Table [T] we give the total number of LBGs found 
in each field in the area covered by both X-ray and optical 
data. 

2.2 X-ray data 

2.2.1 Reduction 

The X-ray data for all the fields were reduced with our 
own pipeline procedure, which uses the Chandra Interactive 
Analysis of Observations (CIAO) software v3.1 and Chandra 
calibration database (CALDB) version 2.27. The reduction 
of the G WS and HDF-N fields is described in lNandra et all 
(|2005bl ) and lLaird etHI (|2005l ) respectively. The same pro- 
cedure was used to reduce the new fields, updating the 
charge transfer inefficiency corrections, applying the most 
recent gain maps, filtering out periods of high background, 
aligning the astrometry between different observations and 
producing images and exposure maps in our various stan- 
dard bands (soft: 0.5-2 keV, fuU: 0.5-7 keV, hard: 2-7 keV, 
uhra-hard: 4-7 keV). 

2.2.2 Point Spread Function 

X-ray source detection requires knowledge of the Point 
Spread Function (PSF) of the Ch andra telescope. W e have 
calculated the PSF using MARX (|Wise et al.ll2003l ). simu- 
lating point sources with monochromatic energy 1 keV at 
a range of positions on the detector. For each position we 
determine the circular aperture which contains a particular 
Enclosed Energy Fraction. The PSF is defined at a position 
on the detector (fixed relative to the mirror). However, our 
deep X-ray data are all the result of a number of obser- 
vations with different orientations and pointing directions, 
merged to create images with the maximum possible expo- 
sure. Thus, to calculate the PSF at positions in our images, 
we determined the EEF radius in each individual exposure. 
We then calculated the circular radius corresponding to the 
average area, weighting by the exposure per pixel in each 
observation. Further details of our PSF calculations will be 
given by Nandra et al. (in preparation) . 

2.2.3 Initial source detection and astrometric correction 

X-ray source detection was initially performed on the soft- 
band image (0.5-2 keV) using the wavdetect tool in CIAO, 
with a low threshold of 10~*. This is used to mask out 
sources when estimating the background in our actual source 
detection procedure (see below). A culled version of the 
WAVDETECT catalogue, containing only higher significance 
detections (src_SIGNIFICANCE> 4), was also used to correct 
for astrometric offsets between the Chandra data and the 
optical imaging. This was done in two steps: first an ini- 
tial cross-correlation of positions was performed using only 



secure, bright optical detections {Ti< 23) from the full SEx- 
TRACTOR catalogue, with a search radius of 5 arcsec; a sec- 
ond matching was then carried out after correcting for any 
offset found, allowing fainter objects {IZ— 23 — 24) to be 
matched within a search radius of 2 arcsec. This allows 
fainter counterparts to also be matched to the X-ray sources, 
increasing the number of matched objects and thus improv- 
ing the accuracy of calculated astrometric offsets. 

2.3 X-ray counts extraction 

Our approach to define a high-redshift sample of AGN is 
based on taking a well-defined, optically-selected sample of 
2 ~ 3 LBGs and determining those which host AGN. Thus, 
rather than searching for optical counterparts to an X-ray 
selected catalogue of objects, within a specified search radius 
we instead extracted soft-band counts at each LBG position 
to determine which LBGs have significant X-ray emission. 
Counts were extracted within circular apertures correspond- 
ing to the 70 per cent EEF for the PSF at that position. 
The background level was estimated locally for each source 
from the number of counts within an annulus of inner radius 
equal to 1.5 times the 95 per cent EEF, and outer radius 100 
pixels greater. All sources detected using wavdetect were 
masked out using 1.5 times their 95 per cent EEF radius 
for the background estimation. The background count rate 
was then rescaled to match the effective exposure and area 
for the source extraction region. Finally we calculated the 
Poisson probability of the source region containing the ob- 
served counts or greater by chance, based on the background 
estimate. An LBG was 'detected' in the X-ray data if this 
probability is less tha n 10~"^. This relative ly high threshold 
(cf. 4 X 10-" used bv lNandra eralll2005b[ ) may be used as 
we have a small number of test positions. Only 0.13 false 
sources in total over all 5 fields are expected due to back- 
ground fiuctuations using this threshold. 

A more problematic contamination, introducing false 
sources into our sample, is due to chance alignment of an 
LBG position with counts from a real X-ray source which 
are not associated with the actual LBG. There are 2 cases in 
which this can happen: i) counts from a nearby bright source 
are scattered into the LBG source region due to the PSF; 
ii) the LBG position falls on top of an X-ray source which is 
not otherwise associated with an optical counterpart. Case 
i) is particularly a problem at high off-axis angles, where the 
PSF is large. To identify these contaminants, we extracted 
counts within the 95 per cent and 50 per cent EEF radii 
at each LBG position. When the LBG position matches the 
X-ray source position, the 95 per cent region should contain 
approximately || — 1.9 times the number of counts in the 
50 per cent region (above background, calculated as above). 
However, if the LBG is close to a bright source, and the de- 
tection is due to scattered counts, this factor will be larger, 
as the 95 per cent region will enclose a larger proportion of 
the bright source PSF distribution. Thus, we identify con- 
taminants due to case i) when the 95 per cent region contains 
more than 4 times the number of net counts as the 50 per 
cent region. 13 significant detections (over all 5 fields) were 
identified as contaminants by this method, and were verified 
as such by visual inspection. 

It is not possible to distinguish between contaminants 
due to case ii) and when the LBG really is an X-ray source. 
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Table 1. Fields used for this work. Column (1): field name; columns (2) and (3): Right Ascension (RA) and Declination (Dec) of the 
centre of the field; column (4): optical filter; column (5): telescope used and date for the optical observations; column (5): seeing FWHM 
of stellar objects in the field (prior to smoothing ) ; colum n (6): total useful optical exposure of data with good seeing; column (7): Galactic 
column density as given bv lDickev fc LockmarJ lll990l ): column (8): X-ray exposure time after good time interval and background flare 
filtering; column (9): survey area with coverage by both X-ray and optical data; column (10) Number of Lyman-break galaxies detected 
in field, in the survey area. 



Field 


RA 


Dec 


Filter 


Telescope and Date 


Seeing 


Optical 


Nh 


X-ray 


Survey 


No. of 


name 












exposure 




exposure 


area 


LBGs 




(J2000) 


(J2000) 






(arcsec) 


(^) 


(lO^Ocm-^) 


(ks) 


(arcmin^) 




(1) 


(2) 


(3) 


(4) 


(5) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


HDF-N 


12:36:55.49 


+62:14:18.28 


K 


Keck I 2003 April 


0.7 


7200 


1.5 


1862.9 


149.1 


292 








G 




0.9 


7560 
















U 


KPNO 2002 April 


1.1 


102600 










GWS 


14:17:43.04 


+52:28:25.20 


K 


KPNO 1996 May/ 


1.2 


8300 


1.2 


190.6 


239.0 


326 








G 


P200 1997 Mar 


1.2 


7200 
















Un 




1.2 


25200 










Lynx 


08:48:55.9 


+44:54:50.0 


n 


WHT 2005 Nov 


1.0 


5500 


2.7 


186.5 


243.4 


223 








G 




1.0 


7000 




















1.1 


12000 










EGSl 


14:22:42.66 


+53:25:24.83 


n 


P200 2005 May 


1.2 


3600 


1.2 


177.8 


358.9 


329 








G 




1.4 


4500 




















1.5 


14400 










LALA 


02:04:44.25 


-05:05:33.83 


n 


WHT 2005 Nov 


1.2 


4500 


2.2 


173.1 


233.3 


144 


CETUS 






G 




1.0 


5000 




















1.3 


22800 











However, we are able to estimate the expected number of 
false sources. To achieve this we generated random LBG 
positions by shifting the original positions by 30-60 arc- 
sec in RA and Dec (this maintains any clustering), and 
repeated our source detection procedure. The process was 
repeated 100 times for each field. The average number of 
false detections in each random sample, after excluding con- 
taminants due to case i) as above, were 1.01, 0.56, 0.46 0.8 
and 0.41 in the HDF-N, GWS, Lynx, EGSl and LALA CE- 
TUS fields respectively. The number varies between fields 
due to a number of factors, such as the depth of the optical 
data and therefore the density of LBGs, depth of the X- 
ray data, area coverage, and PSF effects. However, ~ 70 per 
cent of X-ray sources will have real optical counterparts with 
R < 25, (Georgakakis et al., in preparation), and thus could 
not be identified as LBGs and included in our sample due 
to a chance alignment. The above estimates should therefore 
be corrected by a factor ~ 0.3, and we in fact expect ~ 1.0 
false X-ray detected LBGs in our entire survey. 

The samples of detected objects were then compared 
to previously identified objects in the fields. All but one 
of the X-ray detected LBGs in our sample from the HDF- 
N and GWS have been opti cally identified previously, and 
have spectroscopic redshifts (|Laird et al.ll2006l : ISteidel et al.l 
20031). One objec t has also been optically identified by 
Stern et al.l (|2002l ) in the Lynx field . When no spectroscopic 
confirmation is available, we assign a source a redshift of 
2 = 3, approximately corresponding to the peak of the opti- 
cal selection function (see section [Sjl- Each spectroscopically 
identified source is optically classified as a galaxy (GAL: no 
optical AGN signature), narrow-line AGN (NLAGN: strong 
Lya, significant CIV A1549 emission lines, no emission line 
with FWHM>2000 km s"^), or broad-line QSO (any emis- 
sion line with FWHM > 2000 km s~^). 

To perform the X-ray photometry we use the counts 
within the 95 per cent radius, as the additional counts im- 
prove the statistics. These counts are background subtracted 



and divided by the average value of the exposure map within 
the extraction region, to calculate a count rate corrected for 
instrumental effects (e.g. varying exposure times, mirror vi- 
gnetting, quantum efficiency). This count rate was converted 
to a flux in the same band, and rest-frame 2 — 10 keV X- 
ray luminosities are calculated using the assigned redshift, 
assuming a power-law with F = 1.9 and the appropriate 
Galactic Nh (see Table [T]). 2 — 10 keV rest-frame energies 
are well matched to the observed 0.5 — 2 keV energies at 
z ~ 3, thus large, uncertain fc-corrections are not neces- 
sary. Errors on the fluxes and luminosities are calculated 
from the Poissonian er rors on the to tal observed counts, 
calculated according to iGehreli (|l986l ). There is additional 
uncertainty in the luminosity measurements due to the as- 
sumption of an intrinsically unabsorbed power-law, and the 
assumed redshifts for objects without spectroscopy. These 
effects are small (< 10 per cent) for moderate absorbing 
column densities (Nh < lO'^^cm"^) but could lead to under- 
estimates of the luminosities of more absorbed sources. We 
discuss the possibility of an absorbed population further in 
section [S] We find a total of 22 X-ray detected LBGs. The 
source properties are given in Table (2] 

2.4 High luminosity sample 

Our Lyman-break selected sample of AGN contains very few 
bright objects (Lx > 10*^ erg s~^), due to the comparatively 
small area covered, and we are therefore unable to constrain 
our fit to the XLF at bright luminosities. To remedy this, 
we include bright objects (/x > 5 x 10~^^erg s~^cm~^) from 
the various lar g e area , X-ray selected surveys described by 
iHasinger et all l|2005h . The CDF-N sample is excluded to 
avoid overlap with our own analysis. The X-ray detection 
in these samples was performed in the soft (0.5 — 2 keV) 
band, and is thus consistent with ou r work. Full de t ails o f 
the X-ray data analysis are given bv lHasinger et al] (|2005l ) 
and references therein. 
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Table 2. X -ray detected Lyman-break galaxies. Column (1): Object identification (for GWS and HDF-N adopting designations from 
ISteidel et al ]U003; Laird ct al. 2006; Rcddy ct al. 2006). Column (2): total magnitude in the 7^-band, measured in the Kron aperture 
(see section [2.1.2|l . Column (3) Net X-ray counts above background in the soft (0.5-2 keV) band, within the 70 per cent PSF radius. 
Column (4): soft (0.5-2 keV) X-ray fiux. Column (5): optical classification. QSO - contains broad emission lines (>2000 km s~^), AGN - 
narrow-lined (<2000 km s~^) AGN spectrum, GAL - no signature of an AGN, UID - unidentified spectroscopically. Column (6): redshift 
for objects which have been spectroscopically identified. Column (7): 2-10 keV X-ray luminosity, determined from the soft-band flux 
assuming a power law spectrum with F = 1.9, and the given redshift, or 2 = 3 for UID sources. 



Object ID 


7?. magnitude 


Net counts 


Soft X-ray flux 
(^lU erg s cm ) 


Optical 
type 


Redshift 


X-ray luminosity 
(log_Lx(,^-10 keVj/erg s j 
















HDFN-C14 


24.92 


11.81 


0.55±«:11 


GAL 


2.981" 


42 69+"'io 


HDFN-MD39 


20.50 


1888.95 




QSO 


2.583'' 


44 60+0 01 


HDFN-D21* 


25.37 


3.78 


14+008 


UID 




42.10l°i° 


HDFN-M9* 


24.45 


7.00 


27+0 13 
'-'■^'-0.09 


GAL 


2.975"^ 


42.38t»:i^ 


GWS-M47 


24.93 


8.87 


40+1.08 
■^■^^-0.84 


AGN 


3.026" 


43.50l»:l^ 


GWS-C50 


23.87 


5.70 


2-30to.65 


GAL 


2.910" 




GWS-D54 


22.84 


23.34 


6.83+i?^ 


QSO 


3.199" 


43.86lH«8 


GWS-OMD13 


23.46 


13.87 


4.06tJi« 


QSO 


2.914" 


43.54l0.n 


GWS-MD106 


22.66 


49.83 


15.681?:"! 


QSO 


2.754" 


44 07+0-05 


Lynx- 1 


20.95 


68.84 


26.7612;^ 


QSO 


3.093^= 


44 42+0-04 


EGSl-1 


23.51 


22.87 


1 3 91 +2 '''3 


UID 




44 10+0.08 


EGSl-2 


23.94 


8.85 


9 77+1.12 
^- ' ' -0.82 


UID 




43.40l«:i5 


EGSl-3 


24.09 


18.31 


&.9S+_\l\ 


UID 






EGSl-4 


20.37 


70.94 


on 7r» + 3.62 

00. (U_3 2g 


UID 




44.49tHl 


EGSl-5 


24.15 


7.64 


Q 00 + 1.34 
3.8»_j 02 


UID 




43.55lg:ll 


EGSl-6 


24.36 


10.84 


15.88_3;28 


UID 




44 16+0-10 


EGSl-7 


24.27 


16.21 




UID 




40 Qt: + 0.08 


EGSl-8 


23.79 


20.01 


ii.ool^ig 


UID 




44 00+0-08 


EGSl-9 


23.22 


196.46 


108.09l5?5 


UID 




44 QQ+O-02 


LALA-1 


22.38 


26.50 


13.87l2:« 


UID 




44 in+0-07 

4:4.1U_Q Q7 


LALA-2 


23.75 


7.78 


9 7-1+1.16 
^- '-'--0.84 


UID 




43.39115:11 


LALA-3 


24.27 


4.66 


Qcr + 1.28 
■J-J3-0.95 


UID 




43.49i«:i^ 



Suspected to be starburst galaxies from the X-ray to optical flux ratio (see section [4.41 1. 
Redshift given by Steidel et al. (2003), Laird et al. (2006), " Stern et al. (2002) . 



Spectroscopic identifications in these surveys and at 
these bright fiuxes are highly complete (> 90 per cent). We 
include all objects from th e parent samples with redshifts 
2.5 < 2 < 3.5, whereas the iHasinger et al. I l|2005l ) luminos- 
ity functions were constructed using only unabsorbed 'Type 
1' AGN. In reality, at these bright fluxes and redshifts all 
the AGN in the sample are classifled as Type 1. This may 
be d ue to luminosity d ependence of the absorbed fraction 
(e.g. lSteffen et al.ll2003l ). but the remaining unidentified X- 
ray sources could correspond to a population of high redshift 
Type 2 AGN which will be missing from our bright sample. 
We discuss this issue further in section [S] Our high luminos- 
ity sample consists of 7 objects with Lx = 10^'*'^~*^erg s"^, 
and 5 which are brighter than Lx ~ 10*^ erg s"^. 



3 SELECTION FUNCTIONS 
3.1 X-ray selection function 

The cosmological volume of our survey will depend on the 
depth of our X-ray observations. The sensitivity of Chan- 
dra data varies dramatically over the field of view, due to 
a combination of vignetting and degradation of the PSF. 
To calculate the X-ray selection function we determine the 
number of counts required within a detection cell to sat- 



isfy the detection criterion (false probability > 10~*) for 
each pixel of the X-ray image which is also covered by our 
optical observations. The radius for the detection cell (70 
per cent EEF) and the background are calculated using the 
same method as for the source detection. The count thresh- 
old is then converted to a minimum flux for detection. Thus 
we can calculate the area of our observations sensitive to a 
source above a given X-ray flux. 



For our high luminosity sa mple we adopt the area curve 
given bv lHasinger et ahl (|2005l ) , truncated at our flux limit of 
5 X 10~i^erg s"^. The X-ray detection was performed using 
a variety of methods in the different surveys, and include 
highly significant sources only. However detection is per- 
formed in the same energy range (0.5 — 2 keV), and the vary- 
ing sensitivity limits of these selections are well described 
by the area curve. Thus by incorporating the sensitivity as 
a function of fiux into our maximum-likelihood fitting (sec- 
tion I4.1|l we can correct for the effect of incompleteness in 
the X-ray selection on the survey volume, for the combina- 
tion of our Chandra fields and the various large area surveys 
in a consistent manner. Figure [1] shows the total area curve 
as a function of soft-band X-ray fiux. 
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Figure 1. Total area covered by our survey as a function of lim- 
iting flux. Above 5 X 10~^^ erg s~^cm~^^ ( dotted line) we include 
the surveys used by iHasineer et alj ||2005| ) , which dominate the 
area covered for these brighter fluxes. 



3.2 Optical selection function 

Selecting objects using the Lyman-break technique will al- 
ways miss some fraction of objects at a given magnitude 
due to their true colours not satisfying the selection criteria 
(colour incompleteness) . Photometric errors will also reduce 
the probability of an object being detected at fainter mag- 
nitudes, in addition to scattering objects in or out of the 
selection box. Determining the optical selection function al- 
lows us to correct for these incompletenesses. [Nandra ct all 
ll2005al) used the e ffective survey volumes determined by 
ISteidel et al.l l|l999l ). appropriate for the selection function 
for typical starforming LBGs (i.e. GAL classification above). 
In reality, approximately two thir ds of the X-ray so urces 
have AGN or QSO classifications (|Laird et al.ll200^ . and 
these will have different distributions of intrinsic colours, 
and hence different selection functions. In the present work 
we have calculated the selection functions for each class in- 
dependently. The efficiency of the selection function is also 
highly dep endent on th e dept h and quality of the optical 
data. The ISteidel et al.l l|l999l ) volumes were averaged for 
their set of fields. We determine selection functions appro- 
priate for each individual field in our survey. 

We have calcul ated the selection fu n ctions using Mont e 
Carlo simulations (jSteidel et al.l Il999l : iHunt et al.1 |2004 ). 
The first step was to determine an intrinsic distribution 
of ex pected colours. GAL colours were generated using a 
iBruz ual & Charloi (|2003h model template spectrum for a 
galaxy with continuous star formation and a ge 1 Gyr. This 
templ ate spectrum was reddened using the ICalzetti et al.l 
1I2OOOI ) relation for obscuration by dust, with various ex- 
tinction coefficients drawn from the distribution given in 
lAdelberge r & Steidol (2000) which is based on the observed 
range of LEG spectral shapes. The templates were then red- 
shifted to distribute evenly over the range z = 2 — 4 . Inter- 
galactic absorption was included using the iMadaul (|l995l ) 
extinction curve. Finally the template spectra were multi- 



plied by our filter transmissivities to generate the model 
colour distribution. 

The colours of NLAGN in LBG sampl es are gener- 
ally f ound to be redder than the typical LBG l|Steidel et al.l 
l2002h . We model the expected range of spectral shapes using 
the same approach as for galaxies, but drawing from a Gaus- 
sian distribution of extinction coefficients, roughly based on 
the observed range of UV spectr al slopes for Lyman-break 
selected narrow-line AGN (Steidel et al.ll2002l . 12003, '). Such 
modelling is not a full physical description, being based 
only on a template spectrum for star-forming processes, 
and does not explicitly include any contribution due to the 
AGN. However, it is sufficient to simulate the observed range 
of spectral shapes and redder distribution of broadband 
colours for NLAG N. For QSOs we have followed the work of 
iHunt et al.l l|2004l '). The template spectra were generated by 
varying the continuum slope and Ly g equivalent w i dths o f 
a composite of 59 QSO spectra from lSargent et al] (1 19891), 
based on the Gaussian distributions given in Hunt et all 
(2004). Our model colour distributions are shown in Figure 

El 

Simulations were then performed for each optical clas- 
sification, in which a colour was drawn at random from the 
appropriate distribution, and used to add artificial objects 
to the images for each field. Objects were simulated as point 
sources with seeing appropriate for the field, and 7?,-band 
magnitudes drawn from the range 18.5$;7?.^26. This allows 
photometric errors in the 7^-band to scatter objects into the 
selection range. The objects were recovered using the same 
detection method used for the original data (section 12. ip . 
500 objects were added to each image, and the process was 
repeated 20 times. The selection function, P{z, TZ), was cal- 
culated for a number of magnitude and redshift bins as the 
fraction of input objects which are detected and have colours 
satisfying the LBG criteria. 

The selection functions at bright magnitudes for the 
Lynx field optical data are shown in Figure [3] (a) . Our sur- 
vey varies in efficiency depending on redshift and the optical 
type of the AGN host. The GAL selection function has a 
mean redshift of 2 « 3.0 and samples ~ 65 per cent of the 
total available co-moving volume for a top-hat function for 
z = 2.5 — 3.5. Due to their intrinsically redder colours, the 
NLAGN selection function is less efficient, sampling 50 
per cent of the available volume, and has a lower mean at 
z « 2.8. A strong Lyman-break is not present in the intrin- 
sic spectra of QSOs, but is introduced due to intergalac- 
tic absorption. This results in a higher redshift selection 
function, with a mean redshift z ~ 3.1. The colours are 
also distributed over a smaller range, resulting in a more 
sharply defined selection function, sampling ~ 70 per cent 
of the z = 2.5 — 3.5 volume, but almost 100 per cent for 
z = 2.9 — 3.4. Thus our survey will select a higher pro- 
portion of the overall population of QSOs, biased towards 
slightly higher redshifts, compared to NLAGN. Determining 
the selection functions allows us to account and correct for 
these varying efficiencies. 

At brighter magnitudes (TZ < 22) photometric errors 
are small, and the selection functions vary little between 
fields. At fainter magnitudes, all the selection functions are 
reduced in efficiency and widened, accounting for the larger 
photometric errors increasing scatter, and the increased dif- 
ficulty of detecting faint objects (see Figure [3]d). This is 
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Figure 2. Model colour distributions for different optical classi- 
fications: QSO {triangles), NLAGN {crosses), and GAL {circles). 
Colours are simulated for redshifts 2.0 2: ^ 4.0 (see section[3ll. 
Larger symbols denote objects with 2.7 ^ 2: ^ 3.3. While many of 
the colours falling outside the Lyman-break selection box {solid 
line) correspond to objects at 2 < 2.5 or 2 > 3.5, there is addi- 
tional scatter due to the range of intrinsic spectral shapes. 5000 
model colours are determined for each type (a sample is shown in 
the figure). These colour distributions are used to add simulated 
objects to the optical data, which introduces additional scatter 
due to photometric errors and blending of sources. 



highly dependent on the quality of the optical data in 
each bandpass. For our highest quality data (HDF-N), at 
TZ^i 24.0 — 24.5 the selection function samples ~ 75 per 
cent of the volume sampled at bright magnitudes, and ~ 55 
per cent at 7l~ 25.0 — 25.5. Thus we still sample > 30 per 
cent of the total z — 2.5 — 3.5 volume. In contrast, in the 
EGSl field, with the poorest quality photometry, the effi- 
ciency is 50 per cent that of the bright selection function 
at 7^« 24.0 - 24.5, and ~ 20 per cent at 7^« 24.5 - 25.0, 
for the faintest objects we select in this field. Our incom- 
pleteness corrections are therefore large, but thanks to our 
simulations are well defined over our range of optical mag- 
nitudes. 

Our high luminosity sample is not subjected to a colour 
selection criteria. Instead we rely on spectroscopic follow-up 
of all X-ray sources in the survey fields to have identified a 
complete sample of objects with 2.5 < 2 < 3.5. Thus, we 
assume the 'optical selection' for this sample is a top-hat 
function over 2.5 < z < 3.5 which selects all X-ray sources 
in the redshift range. This is a reasonable assumption given 
the very high completeness of the samples at bright X-ray 
fluxes. This can then incorporated into our analysis in an 
analogous manner to the incomplete selection functions for 
the Lyman-break sample. 
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Figure 3. Optical selection functions using the Lyman-break 
technique, shown here for the Lynx field optical data. The redder 
colours of NLAGN {dashed line) result in a less efficient selection 
function shifted to lower redshifts than for galaxies {solid line). 
The lack of an intrinsic Lyman-break in the spectrum of QSOs 
moves the selection function to higher redshifts {dotted line). The 
colours of QSOs are distributed over a smaller range than for 
galaxies or NLAGN, leading to a more sharply defined selection 
function. For bright objects (a), the main source of incompleteness 
is the colour selection. At fainter magnitudes (b), the selection 
functions are less efficient due to detection incompleteness and 
larger photometric errors, and is dependent on the depth and 
quality of the optical data for a field. 



4 LUMINOSITY FUNCTION 

A great variety of methods have been used in the past to cal- 
culate luminosity functions. One of t he simplest an d most 
widely used is the 1/Vmax method l|Schmidtl Il968l ). Each 
object is assigned a maximal survey volume corresponding 
to the redshift interval over which the source could be in- 
cluded in the sample, given a redshift selection range and the 
flux limits of the survey. Usually binned luminosity functions 
are constructed by summing the contribution of all objects 
within a luminosity bin: 



E 

AlogLx 



However, significant errors are introduced by this method 
when the objects are close to the flux limits l|Page &: Carreral 
l2000h . and it is poorly suited for low numbers of sources. Our 
sample is small, and the flux limits of our survey vary widely 
between different fields, as well as over a single field of view. 
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We therefore take a different approach, using a maximum 
likehhood fitting method. No binning is required, so infor- 
mation on individual sources is retained. We are also able to 
include completeness corrections, correcting for the varying 
sensitivities of both the X-ray and optical data. In section 
I4.2l we calculate improved binned estimat es of the luminosity 
function using the i Vobs/A/mdi method (jMivaii et al.l I2OO1I : 
iHasinger et al.ir2005h for display purposes and comparison 
with other work. 



4.1 Maximum likelihood fitting 

To perform a maxi mum likelihood fit we use the STY 
l|Sandage et al.ll 19791 ) estimator, which requires an assumed 
functional form for the luminosity function. We adopt a bro- 
ken power-law for the differential luminosity function, which 
has been found to descr ibe the shape well at lower redshifts 
(e.g. iBareer et ai]l2005l ). 



d$(Lx; 
d log L> 



+ 



Lx 



(3) 



where Lx is the 2-10 keV (rest-frame) X-ray luminosity. The 
probability of observing a source of a particular luminosity, 
Lx, is then given by 



Pi = 



I42 I2.5 4>{Lx) A{Lx, z) C{Lx, z) ^ Az d log Lx 



(4) 



The denominator is the expected number of sources at z ^ 3, 
found by integrating </>(Lx) over the survey volume. A(Lx, z) 
is the total area from all our fields sensitive to the fiux from a 
source of luminosity Lx at redshift z. The luminosity to flux 
conversion assumes an unabsorbed power-law with F = 1.9 
and the appropriate Galactic Nh for each field, and thus is 
consistent with our original calculation of luminosities from 
observed flux, ^dz is the differential co-moving volume el- 
ement per unit area. 

Two further corrections for incompleteness effects are 
included in our calculations. In Equation |4] the factor 
C(Lx, z) reduces the expected number of sources, account- 
ing for objects fainter than the optical magnitude limit. Fig- 
ure |4l which shows the optical TZ magnitude vs. soft band 
flux illustrates the need for this correction. Lines of constant 
X-ray to optic al flux ratios are plotted, ac cording to the re- 
lation given in iHornschemeier et al.l ll200ll). converted for 71 
magnitudes using the relation in lSteidel fc HamiltonI (|l993l ) 
(with no optical colour correction term). 



(5) 



For AGN, this ratio is distributed over the range —1.0 < 
log (/x/ /opt) < -1-1.0 over s evera l orders of mag r iitude 
in flux (e.g. ISchmidt et al.1 1 19981 : lAkivama etaU I2OO0I : 
[Alexander et al.ll200ll ). Our sample of X-ray detected LBGs 
fall within this range. However, at faint X-ray fluxes (/x < 
10~^^erg s~^cm~^) we would expect some fraction of the 
sources to have optical magnitudes TZ> 25.5, and so would 
not be found using our Lyman-break selection. The C{Lx, z) 
factor in equation|4]attempts to correct for this, by reducing 
the expected number of sources in proportion to the fraction 
of the —1.0 < log(/x//opt) < +1-0 falling below our optical 
magnitude limit, and assuming an even distribution over this 



range. Some small f raction of extreme AGN may reside out- 
side this range (e.g. Hornschemeier et al.ll20ol] : iBarger et al.l 
I2OO3I : iLehmer et al.ll2005l ). but the true distribution is un- 



known, necessitating this simplifying assumption. Thus we 
adopt, 

CiLx,z) = ^ + 1(1.0 + log /x(Lx,^) + 5.4) (6) 

This correction is directly related to the X-ray flux, and 
thus is defined for every value of Lx and z. For fiuxes ~ 
10~^^erg s~^cm~^, corresponding to the median flux in our 
sample, we expect 20 per cent of the sources to have TZ > 
25.5, so this correction is comparatively small, but rises for 
the fainter fluxes probed by our survey. At the limit of our 
~ 200 ks fields, /x ~ 2 x 10~^'^erg s~^cm~^, the correction 
is a factor ~ 2. At the extremely faint X-ray fiuxes probed in 
the HDF-N we actually expect all the optical counterparts 
of AGN to be fainter than our limit. Thus, we suspect the 2 
faintest objects in the HDF-N have a starburst origin, which 
we discuss further below, and exclude them from our fitting 
of the XLF of AGN. 

Our second correction accounts for incompleteness 
within the Lyman-break samples (i.e. at 7l< 25.5). The 
optical selection functions calculated in section |31 even at 
bright magnitudes, are considerably less efficient than a 
z = 2.5 — 3.5 top-hat function, implicitly assumed in the 
integration limits in Equation!?] However, this correction de- 
pends on the optical classification and 7?.-magnitude, which 
is specific to each detected source. Thus, we correct for this 
effect by introducing an incomplet e ness weighting, Wi, for 
each object, following IZucca et al. I (fl994i ). Wi is inversely 
proportional to the probability of object i being included 
in our sample, and thus can be calculated from the ratio 
of the total survey volume (given the X-ray luminosity of 
the source, Li), and the effective volume sampled by the 
selection function. 



Kot 



2.5 ML^,z)—dz 



j:^.MsIML^,z)Pf{n,,z)^dz 



(7) 



where Pf{lZi,z) is the selection function for a field at the 
objects optical magnitude, TZi, for the optical classification 
type, as determined in section|3l and Af{Li, z) is the area of 
that field sensitive to the source luminosity, Li, at z. We sum 
the effective volumes from each of our fields. For spectro- 
scopically unidentified objects we must use an average of the 
three selection functions, weighted acc ording to the f r action 
of each spectroscopic type expected. 'St eidel et al.l (|2003l ) 
find roughly equal numbers of LBGs classified as NLAGN or 
QSO in the ir very large samp le of spectroscopically identi- 
fied objects. iLaird et al.l (|2006l ) also found roughly one third 
of X-ray detected LBGs exhibited no optical signature of 
an AGN. Thus we take the average of the three selection 
functions, with each type weighted equally, for unidentified 
objects. For our high luminosity sample Pf{TZi,z) is a top- 
hat function, and thus each object is assigned a weighting of 
Wi = 1. Figure |S] demonstrates the effects of our complete- 
ness corrections on the determination of the faint-end slope 
of the XLF. 

The best fit values for 71, 72 and L« in equation [3] can 
then be found by minimizing the log-likelihood function. 



S = -2\nC = -2^ 



<w > 



In Pi 



(8) 
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Table 3. Best fit parameters for the XLF at z < 
the maximum UkeUhood fitting (section 14. 



' 3 derived from 



Parameter 



Value 



71 
72 

log L, / erg s"-' 

-ftTnorm / MpC'S 



54+0-20 

r, oq+0.40 
^■■'•^-0.30 
44.68i°:f 

4.58ti;oI X 10-6 



Introducing the weights artificially reduces the size of errors 
estimated from the likelihood fitting, and is th erefore bal- 
anced by the average weight jllbert et al.ll2005l ). This does 
not effect the minimization of S to determine the best fit 
parameters. 

The overall normalisation of the luminosity function is 
not determined from the minimization of equation |8] and so 
is calculated using the best fit model, 0mdi(ix), to match 
the observed number of sources, 



I /2^5^ <?^mdl(ix) A(Lx,z) C{Lx,z)^ dz dlogLx 

(9) 

where Wi gives the effective total number of sources, cor- 
rected for incompleteness in the optical selection. 

We estimate la errors on the model parameters from 
the maximum projection of the AS* =1.0 likelihood surface 
for each parameter. The error in A'norm is estimated from 
the Poissonian error given the total number of sources. The 
best fit parameters are given in Table [3] 



4.2 Binned estimates 

Binned estimates of the luminosity function a re calculated 
using the Nobs/N^^i method described by iMivaii et al.l 
l|200lf ). The value within a bin is given by: 



"^bin — <^mdl(-£'bin) 



iVol 



(10) 



where ^mdi is the best fit model from the maximum likeli- 
hood fitting, evaluated at the bin centre, Lbin- A^'obs corre- 
sponds to the observed number of objects within the bin, 
but for our analysis this should be corrected for the Lyman- 
break selection incompleteness. Thus, 

Nohs = 

A log Lx 

N^di is the expected number of objects within the bin based 
on the best fit model. 

^3.5 
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Figure 4. X-ray flux (0.5-2 keV) vs. TZ magnitude for the sample 
of X-ray detected LBGs. Squares: QSO optical classification; cir- 
cles: NLAGN; triangles: GAL; inverted triangles: no optical iden- 
tification. The majority of our sources are found to be distributed 
between —1.0 < log(/x//opt) < +1.0, the expected range asso- 
ciated with AGN. The Lyman-break selection limits the faintest 
7^-magnitude at which we can detect a source to TZ< 25.5 (dotted 
line), necessitating the incompleteness correction given by equa- 
tion |6] The two faintest objects in our sample (s tars) are sus- 
pected to be star-burst galaxies l lLaird et al.ll200g . section |4j4jl . 
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This estimator accounts for the varying sensitivity limits 
within a bin, which are known to re sult in biased est i mates . 
It offers an improvement over the IPaee fc Carreral (|2000l ) 
(peat method as it allows the luminosity function to vary 
within a bin. Errors in the binned estimates are calculated 
from the Poisson uncertainties for the number of objects in 
each bin. The binned results, along with the best fit, are 
plotted in Figure [S] 



Figure 5. X-ray luminosity function of AGN at z ~ 3. Black line: 
our maximum likelihood fit; black circles: binned estimates us- 
ing the Afjj^js /A''nj(ji method; open orange circle: l/Vmax estimate 
for the lowest luminosity bin, containin g the s uspected star-burst 
galaxies; red dashed line: iBarger et al. I 1I2OO5I) best fit for a PLE 
model at z 1.2, evalua t ed at z = 1; green triangles: reproduc- 
tion of the iBarger et al. I l l2005h work/method for z = 2.5 - 3.5, 
using spectroscopic and photometric redshifts. 
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Figure 6. The effect of the completeness corrections on the XLF. 
Solid black line/ circles: our best fit XLF and binned estimates, 
with all completeness corrections applied; dashed blue line/open 
circles: without accounting for objects fainter than 1Z = 25.5 
based on the /x//opt distribution; dotted red line/open squares: 
no correction for optical incompleteness (a top-hat function over 
z = 2.5 — 3.5 is assumed for all optical magnitudes, rather than 
the Lyman-break selection functions, and there is no correction 
for objects fainter than = 25.5). Poissonian error bars are only 
shown for the best fit, but will be identical for the other points 
at the same luminosity. Correcting for the incompleteness of the 
Lyman-break selection increases our determination of the XLF 
at the faint-end, but does not steepen the slope. Accounting for 
objects fainter than TZ = 25.5 results in a st eeper faint-end slope . 
Open green triangles: reproduction of the iBarger et al.l l|2005l ) 
work/method for z = 2.5 — 3.5, using only (incomplete) spectro- 
scopic identifications. 

4.3 Evolution of the X-ray luminosity function 

Our z = 3 XLF shows excellent agreement with the double 
power-law model above a luminosity of log Lx = 42.5. Be- 
low this we see an upturn, which we comment further upon 
below. For comparison wi th lower redsh i fts, w e plot the best 
fit XLF on Figure [5] from iBarger et al] (|2005l ) , evaluated at 
z = 1 for a PLE model fit to objects with z < 1.2. It can 
be seen that both the bright and faint slopes a,t z — 3 are 
similar to those at z — 1, but the XLFs are clearly differ- 
ent in that the z — 3 data points at low luminosity fall 
below those at z = 1. In addition, L, at z = 3 appears 
to be systematically higher. To make a formal comparison, 
we have determined two-dimensional confidence regions for 
71 and which we show in Figure [71 compared to the 
IBarger et all (|2005l ) parameters for the PLE model. Due to 
our small sample size, the errors in our best fit parameters 
are large and dominate over the errors at z = 1. None the 
less, our determination of 71 is in good agreement with the 
Barger value at lower redshifts. We thus find no evidence 
for a fiattening of the faint-end slope of the XLF between 
redshift 2 = 1 and z ^ 3, although we are unable to exclude 
the possibility of some amount of flattening. We do, how- 
ever, flnd weak evidence for continued positive evolution of 
I/t , but the significance of this result is not high. We find no 
evidence for evolution i n the bright-end, although we note 
that the IBarger et al.l l|2005l ) 2 = 1 XLF and our 2 ~ 3 
results are poorly constrained at Lx ^ lO^^'^erg s~^. 



Figure 7. Confidence interval for L, and 71 for ^ ~ 3 maximum 
likelihood fitting {solid: 68.3 per cent, dashed: 95 per cent confi- 
dence le vel, error bars i ndica te 1 dimensional 1 a errors). We also 
show the I Barger et all ||2005|) PLE model, evaluated at 2 = 1 (tri- 
angle) . We find no evidence for a flattening of the faint-end slope 
(although we cannot rule out the possibility), and weak evidence 
of continued positive evolution of L,. 

Integrating the luminosity functions, we find a total 
number density of AGN (with logLx > 42) of 1.0j:o;4 erg 
s~^ Mpc~'^, a decline by a factor ~ 3.6 from 2 = 1 to 2 = 3. 
Because of the increase in L* , the decline in the inferred lu- 
minosity density is much smaller. We find a total 2-10 keV 
luminosity density for black hole accretion of 2.21q 4 x lO'^^ 
erg s~^ Mpc~^, about a factor ~ 2 smaller than at 2 = 1. In 
Figure |8] we show results for the space density derived from 
our best fit XLF at 2 ~ 3 in two luminosity ranges. There is 
a clear decline in the space density of these lower luminosity 
AGN ab ove redshift 2^ 1, and our results are consistent 
with the lUeda et"aD (|2003l ) LDDE model. 



4.4 Pure starburst contribution 

The two faintest X-ray sources in our sample, both detected 
in the HDF-N, have an X-ray to optical flux ratio which 
places them outside the expected range for AGN. We there- 
fore suspect that the X-ray emission from these objects is 
due to vigourous starbursts taking place, rather than AGN 
activity. It is not unexpected that the only detections at 
these faint fluxes are starbursts, rather than AGN. The TZ- 
band limit for initial pre-selection with the Lyman-break 
technique means we do probe deep enough with the opti- 
cal data to sample the AGN /x//opt range (see Figure 3)). 
One of suspected starbursts (HDFN-M9) was previously re- 
ported by Laird et al. (2006). This object has been spectro- 
scopically identified, and is classified as GAL due to the lack 
of any signature of an AGN in the optical spectrum, which 
is consistent with a starburst origin for the X-ray emission. 

These objects are contaminants of the AGN sample, 
and were thus excluded in our fitting of the XLF. It is how- 
ever informative to calculate the differential space density 
in our faintest bin, including these 2 detected objects. The 
/x//opt completeness correction (equation |6]l is not valid for 
these objects since they are not AGN, and we cannot use the 
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Figure 8. Space density of AGN as a function of redshift, 
for Lx(2 ~ 10 keV) = lO'^^-'^'^ergs-^ {solid line) and Lx(2 - 
10 keV) = 10^*~^^ erg s~^ (das/ied line). Lines show the best fit 
Lumin osity Dependent Density Evolution model of lUeda et al.l 
l|2003l V in which the number density of AGN peaks at z ~ 1, 
with the peak moving to lo wer redshifts for lower luminosities. 
We also show the res ults fromlBar ^er et al. (2005) at low redshift 
{open squares), and iBareer fcCowiel 1I2OO5I) for 2 = 2 — 3, us- 
ing IR and optical photometric redshifts to improve completeness 
{solid squares). Our results { circles) using Lyman-b reak selection 
are in good agreement with iBarger &: Cowig (|2005l ) , and consis- 
tent with the decline in number density of lower luminosity AGN 
above 2 ~ 1 predicted by LDDE. 



Nobs/N^di method to calculate the binned estimate as this 
is dependent on the fit to the AGN XLF. We therefore use a 
simpler 1/Vcs estimate, where Vcs is calculated as in equa- 
tion [71 and we sum the contribution of the 2 objects for a 
0.5 dex luminosity bin. This result is shown in orange in 
Figure O and is considerably higher than the AGN XLF. 
Comparing to the predicted number of AGN detectable in 
our survey, based on our best fit XLF, we estimate the sig- 
nificance of this excess at > 2a. 

This upturn in the XLF at the faintest luminosities has 
not previously been seen at 2; ~ 3, as preceding studies 
have not probed to sufficiently faint fiuxes. Neither of our 
suspected sta rburst galaxies were e ven reported as X-ray 
sources m the lAlexander et al.l (120031 ) catalogue of the CDF- 
N. Given we only detect 2 objects, there is considerable un- 
certainty in our result. The formal significance of the excess 
contribution is low, and we are limited to a very small area 
with sufficiently deep X-ray data. We may also be particu- 
larly susceptible to Eddington bias so close to our detection 
limit. However, detection of an additional suspect ed star- 
burst galaxy was also reported bv lLaird" et al.l(|2006l ). which 
now falls just below the detection limit when using our up- 
dated PSFs and direct extraction of the X-ray counts. We 
thus conclude that at high redshift, such starburst galaxies 
are likely to be a significant contaminant of the XLF of AGN 
at t he very fa intest luminosities. 

iLaird et'al. (2006) performed a stacking analysis of the 
X-ray counts from the remaining LBGs in the HDF-N with- 
out direct X-ray detections, which allows us to probe to 
yet fainter X-ray luminosities, below our detection limit. 
We have derived the space density of these objects, using 
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Figure 9. Extrapolation of the faint-end of the integrated lumi- 
nosity function at 2; ~ 3 for comparison with the contribution of 
star-forming galaxies. Both the 1/Vmax estimate for the number 
density of the 2 suspected star-burst galaxies {open circle) and the 
estimated number density of star-forming galaxies { square, from 
the stacking analysis of X-ray undetected LBGs by iLaird et al.l 
fall above our best fit luminosity function. 

the 1/Vmax method and our GAL selection function at the 
7?,-magnitude of each object. In Figure|9]we extrapolate our 
best fit XLF to lower luminosities, and compare to the space 
density of the stacked LBGs at their average X-ray luminos- 
ity, and our 2 starburst detections. The stacking result also 
falls above the predicted contribution from AGN. Combined 
with our direct detections, this provides further evidence for 
a large population of bright star-forming galaxies at high 
redshift, which dominate in number density over AGN at 
these luminosities. This contribution has not been found 
in lower redshift {z < L2) studies (e.g. iBarger et al]|2005l ) 
of the XLF, indicating continued evolution of the bright- 
est star-forming galaxies above z = 1, consistent with find- 
ings from UV /opt ical studies (e.g. ISchiminovich et al.ll2005l : 
iReddv et al]|2008l ). 



5 DISCUSSION 

5.1 Evolution of AGN at high redshift 

Our results show that the XLF of AGN continues to evolve 
above redshifts z ^ 1. This evolution is characterised by 
a reduction in the number density at faint luminosities 
(Lx < I/*), although the actual faint-end slope may not flat- 
ten. Such evolution cannot be explained as pure luminosity 
evolution. The faint-end evolution appears more consistent 
with a reduction in the density. However, there is clearly a 
luminosity dependence as this evolution is not seen in the 
bright end. Our result is thus consistent with some form of 
Luminosity Dependent Density Evolution. However, most of 
the specific forms of LDDE in the published literature in- 
volve a flattening of the faint end slo pe above z ^ 1 (e.g. 
iHasinger et al. I I2OO5I : lUeda et al.|[20o3 ). We find no such ev- 
idence, and it is possible that the inference of a flattening 
seen in previous work is due to incompleteness in the redshift 
identifications of X-ray sources. While our Lyman-break se- 
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lection also suffers from incompleteness, we have corrected 
for such effects in our calculations. 

The initial results for the space density of moderate lu- 
minosity A GN (Lx = 10' * '^"'*'*'^ e rg s~^) at z ~ 3 using our 
method, bv lNandra et al.1 l|2005al ) , showed a markedly higher 
value than previous estimates, particularl y those relying pre- 
dominantly on sp ectroscopic redshifts (jCowie et al.l l2003l : 
iBarger et al."2005 !). Ou r new results are consistent with both 
those of iNandra et al.1 l|2005al ) an d the revised space den - 
sity us ing photometr i c redsh ifts bv lBarger fc Cowi3 l|2005l ). 
While iNandra et al.l (|2005al ) found little evidence for a de- 
cline in AGN numbers from z = 1 to z; = 3, the improved 
analysis and larger sample in our current work strongly in- 
dicates that there is such a decl ine, in support of the conclu- 
sions of iBarger fc Cowid l|2005l ). Indeed, our results should 
be more robust than either of these previous studies, due 
both to the improvements in the combined analysis of X-ray 
sensitivity limits and the optical selection functions, and the 
much larger survey area and source numbers at faint fluxes. 

We can also comp are our results with p revious estimates 
of the XLF at z ~ 3. IBarger et all HqoI) presented high- 
redshift XLFs in two bins, z = 1.5 — 3 and 2 = 3 — 5, which 
are not well matched to the redshift range sampled by our 
Lyman-brea k selection. We have therefore attempted to re- 
produce the IBarger et al.l l|2005l 'l XLF using the same sam- 
ples, but for z = 2.5 — 3.5. We select soft band detected 
sources wit h spectroscopic or photome tric redshifts from 
the CDF-N (|Barger et al.ll2003l ). CDF-S (| Zheng et al.ll2004h . 
and Chandra Larg e Area Synoptic X-ray Survey (CLASXS, 
ISteffen et al.ll2004l). The area covered as a function of flux is 
taken from^Bar ger et al] (|2005l ). The 1/Knax method is used 
to calculate the binned XLF, thus. 
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where the summation is over all A^'bin sources in the lumi- 
nosity range for each bin. 

Initially we only include sources with spectroscopic red- 
shifts. The results are shown by the green triangles in Figure 
\6\ While at luminosities > lO^'^'^erg s~^ there is good agree- 
ment with our best determination of the z ~ 3 XLF, at lower 
luminosities there is only agreement with our results if we 
do not apply completeness corrections. Our corrected result 
is considerably higher. Including objects with photometric 
redshifts increases the numbers at lower luminosities, and 
there is much better agreement with our XLF (see Figure 
O, although our result is still systematically higher. Clearly 
using incomplete samples of low luminosity X-ray sources at 
these redshifts, without applying completeness corrections, 
can incorrectly indicate a flattened faint-end slope. Photo- 
metric redshifts are necessary to increase the completeness, 
but at the faintest luminosities there may be residual in- 
completeness which is difficult to correct for. There remain 
a number of other uncertainties with photometric redshift 
determinations, including the errors in the photo-z's, pos- 
sibility of catastrophic failures, and the high probability of 
chance superpositions of optical counterparts. 

We also note the clear co-dependence between determi- 
nations of 7i and L, in Figure [T] Accurate measurements of 
the faint-end slope are therefore essential to determine L* 
accurately, which is often used as a tracer of the more funda- 



ment al underlying black-h ole /host galaxy mass distribution 
(e.g. iHopkins et al.ll2006l ). Thus, careful completeness cor- 
rections for the detection of faint AGN, to ensure that the 
faint-end slope is not artificially flattened, are essential when 
using the luminosity function to investigate AGN evolution. 

While we find marginal evidence for positive evo- 
lution of L,, we find no evidence for evolution of the 
XLF at brighter luminosities (Lx ^ lO^^erg s~^), which 
has previously been suggested by a number of authors 
fe.g.lUeda et al.ll2003l : iHasinger et"al]|2005l : lEa Franca et al.l 
l2005h . This may be purely due to uncertainty in our bright- 
end determination, where we have small numbers of sources. 
A high luminosity, high redshift population may be missed 
due to remaining incompletness in the spectroscopic follow- 
up of our high luminosity sample, although this will also 
affect previous studies. There is evidence that the number 
density of the brightest AGN peaks a round z » 2, but de- 
clines again at highe r redshi fts (e.g. IHasinger et al ] I2OO5I : 
ISilverman et~aLl |2005| . l2007a|)^ in a similar manner to op- 
tically selected QSOs ^e.g. IWolf et~aLl l2003l : ICroom et all 
I2OO4I : iBongiorno et al ] l2007l ). which could explain the ap- 
parent lack of evolution when comparing z = 1 and z — 3. 
The lack of observed evolution may indicate that we miss a 
population of luminous, absorbed AGN (i.e. Type 2 QSOs) 
which are detected in the hard band at z = 3. The lack 
of AGN in our high luminosity sample which exhibit opti- 
cal signatures of absorption supports this conclusion. Such 
a population may not be identifiable at z = 1 as the hard 
band probes lower rest-frame energies. 



5.2 Limitations of the LBG selection method 

We believe our results using the Lyman-break method to 
select, and correct for incompleteness in high redshift sam- 
ples of X-ray detected AGN, are the most robust and ac- 
curate yet, but there are limitations to our approach. The 
necessary completeness corrections are well defined (an ad- 
vantage over direct identification methods), but they can be 
quite large, particularly for the faintest X-ray luminosities. 
At the Lx ~ 10'*^~'*^'^erg s~^ our typical overall correction 
is a factor ~ 4, accounting for both the incompleteness in the 
Lyman-break selection function for the typical 7^-magntude, 
and the expected fraction of sources fainter than our optical 
limit. This incompleteness means we require a larger area 
coverage to obtain comparable sample sizes to directly iden- 
tified samples of high redshift X-ray sources, although the 
incompleteness of such samples will also be high, and more 
difficult to correct. 

In section 13.21 we derived selection functions for differ- 
ent optical classifications, and noted the differences in red- 
shift and overall efficiency. However, we lack spectroscopic 
confirmation for more than half our sample, and thus use 
an average of our 3 selection functions. Spectroscopy would 
also allow us to determine the true redshifts of our uniden- 
tified sources, which are currently assigned z = 3. However, 
while spectroscopy is clearly desirable, the errors in our re- 
sults are dominated by the Poissonian uncertainties due to 
the small number of sources, rather than redshift or optical 
classification. 

A number of ot her authors (e.g. lUeda et al. I I2OO3I : 
iLa Franca et al.ll2005h have attempted to correct for the ef- 
fect of intrinsic absorption, due to obscuration within the 
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AGN host, on the X-ray luminosity. We have chosen not to 
do this for a number of reasons. Our work is focussed on 
2^3, thus the observed soft (0.5 — 2 keV) emission corre- 
sponds to a comparatively hard rest-frame energy (~ 2 — 8 
keV), which should be reasonably unaffected by moderate 
absorbing column densities (Nh < lO^'^cm"^), and so correc- 
tions should be small. Many of our sources contain very few 
counts, thus there is little spectral information and absorp- 
tion corrections are uncertain. Additionally, if an absorp- 
tion correction is applied this removes the simple relation 
between observed flux and luminosity, which would compli- 
cate our analysis in section 371] Performing X-ray detection 
and photometry at harder observed energies, but with the 
same analysis, would be possible to reduce absorption ef- 
fects, but would substantially reduce the sensitivity to faint 
fluxes. Indeed, only 11 of our soft X-ray detected LBGs are 
also detected in the hard band to the same significance level. 
None the less, it remains a limitation of our study that we 
are unable to correct for absorption effects without running 
the risk of introducing greater error. 

5.3 Census of high redshift AGN 

The comparison between our results and previous estimates 
of the XLF at high redshift provides an indication of the 
nature of high-redshift AGN, and their host galaxies. The 
faint-end of our Lyman-break selected XLF is systematically 
higher than determined from (incomplete) foUowup of purely 
X-ray selected samples. Thus the Lyman-break method se- 
lects a sample of objects which accounts for the entire known 
population of X-ray selected, moderate luminosity AGN at 
2: = 3. To be selected, these AGN must be bright in the rest- 
frame ultraviolet, and have similar colours to star- forming 
LBGs. This is hardly surprising for bright QSOs, whose UV 
emission will be dominated by the AGN. For fainter objects, 
even those with spectra, it is difficult to determine the rela- 
tive contribution to the UV flux from the AGN accretion ac- 
tivity, and the s tar-formation proces ses taking place within 
the host galaxy. ISteidel et al] (|2002l ) examined a sample of 
optically identified AGN in LEG samples, and concluded 
that the host galaxies of observed AGN could be the equiv- 
alent of LBGs without active nuclear activity. Using X-ray 
emission to identify AGN in the LBG samples finds an ad- 
ditional population which do not exhibit optical AGN sig- 
natures, and supports the conclusion that z ~ 3 AGN are 
hosted by galaxies with similar properties to the LBGs. 

This is of great interest, as the majority of moder- 
ate luminosity A GN at lower redsh i fts are found in re d 
host galaxies fe.g. iBarger et al.ll2003l : iNandra et al]l2007l l. 
iNandra et all (|2007l7 presented a colour-magnitude diagram 
for a sample of X-ray selected z ^ 1 AGN, which were found 
to lie mainly in moderate luminosity galaxies (—20.5 < 
Mb < 22.5) on the red sequence, at the top of the blue 
cloud and in the valley in between. 90 per cent of the X-ray 
sources have rest-frame U — B > 0.6. By contrast, the ab- 
solute magnitude limit of the LBG samples is Afs ~ 22 at 
z = 3, although our /x//opt completeness correction does 
attempt to account for AGN residing in fainter hosts. 

We are unable to directly compare the z ~ 1 colours 
to our sample, as the rest-frame U and B bands correspond 
to much redder wavel engths at z ~ 3 than covered by our 
photometry. However, IShaplev et all (|2005l ) have performed 



SED fitting for a sample of rest-frame UV colour selected 
galaxies, utilising ground based photometry from 0.35 to 
2.15 (xm, complemented by Spitzer IRAC data. We have de- 
rived rest-frame U — B colours for each of these galaxies, 
based on the best fit star formation histories. The mean 
{U — B)ieat colour is 0.52, with standard deviation 0.23. All 
of the X-ray selected 2: ~ 1 AGN (excluding QSOs) are 
redder than this mean, and 52 per cent lie more than 2 
standard deviations away, and are thus redder than over 98 
per cent of the rest-frame UV selected galaxies. Assuming 
the host galaxies of our z ~ 3 X-ray detected LB G sam- 
ple have similar colours to the lShaplev et all lj2005h sample, 
this tentatively suggests a significant evolution of the host 
galaxy colours from the blue cloud at 2 = 3 to predomi- 
nantly the red sequen ce and green valley ai z — 1 (see also 
[Silverman et al.ll2008l ). 

However, this result does not exclude the possibility of a 
population of moderate luminosity AGN residing in red host 
galaxies at 2 = 3, but they must be under-represented in 
both Lyman-break and directly identified samples. Alterna- 
tive colour selection techniques have identified populations 
of red ga laxies at high redsh ifts (e. g. BzK passively e volving 
galaxies, iDaddi et al1l2004l : DRGs, iFranx et al]|2003l ') which 
have little overlap with rest-frame UV selected samples 
iReddv et al.ll2005l ). Such galaxies could host faint, mod- 
erately obscured AGN. Our correction for objects fainter 
than the 7?.-band limit may partially account for such a 
population, which will be faint in the rest-frame UV, but 
it is not possible to determine if they are adequately rep- 
resented. Thus, if there is a significant population of AGN 
in red galaxies it will under-represented in our work. Direct 
spectroscopic identification or photometric redshift deter- 
minations will also be difficult due to the faintness of such 
galaxies, and thus they could also correspond to the uniden- 
tified fraction of purely X-ray selected samples. 

As discussed in section 15.11 hard-band selected samples 
of high redshift AGN may include an absorbed, high lumi- 
nosity population which is missing from our soft-band se- 
lected work. Absorption will have an even more severe effect 
on the faint-end of the XLF, as moderately absorbed sources 
will not be identified in the soft band, and the hard band will 
not probe deep enough. Thus a significant population of low 
luminosi ty absorbed AGN m ay r emain unidentified a t high 
redshift. lOaddi et all (|2007l ) and iGeorgakakis et all (|2008[ l 
have reported populations of faint, hard sources based on 
stacking analyses, but the key issue of whether they are 
Compton-thick is still unresolved. 

5.4 Future prospects 

While our survey covers the largest area of deep (> 200 
ks) X-ray data, and probes to the faintest limits, our total 
sample is still small, and thus uncertainties remain. Fur- 
ther work, increasing the sample size, is required to im- 
prove constraints on the parameters describing the XLF at 
2 ~ 3. The Lyman-break technique is extremely efficient for 
identifying AGN in large area X-ray surveys, such as those 
covered by recent ~ 200 ks Chandra surveys (e.g. E-CDFS, 
iLehmer et ahlbOOSl. AEGIS-X , PI: Nandra; C-COSMOS, PI: 
Elvis). IPrancke et al.l ()2008l 'l have already applied similar 
techniques in the E-CDFS to identify moderate luminosity 
AGN at 2 ~ 3 and perform a clustering analysis. Combined 
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with well-defined completeness corrections these deep, large 
area surveys should vastly improve the constraints on the 
faint-end (Lx « 10''^"'''*-^erg s"^) of the XLF. Improving 
constraints at > L« requires even larger area coverage of 
more moderate depth, and thus the efficiency of our Lyman- 
break method may also be beneficial for such studies. A 
larger area of data of depth > 2 Ms or greater is required to 
probe to Lx < 10^'^'^erg s~^ and confirm our X-ray detec- 
tion of a starburst population. However, forthcoming deeper 
X-ray observations of the AEGIS field (A09, PI: Nandra), 
taking 3 of the current fields to a depth of 800 ks, will vastly 
increase the area sensitive to Lx > IC^^'^erg s~^, without 
the issue of starburst contamination, and provide essential 
data to improve constraints at the faintest fluxes, currently 
only accessible by the Chandra Deep Fields. However, to 
really constrain the high-redshift evolution of the XLF re- 
quires large samples, selected over a range of redshifts, and 
probing to low luminosities, requiring the combination of 
all surveys over the full range of depths. However, to fully 
utilise the available depth will require carefully derived com- 
pleteness corrections, which must therefore be an essential 
feature of any future studies. 



6 SUMMARY 

We have presented a new determination of the XLF of AGN 
at 2 = 3, by combining Lyman-break selection in the optical 
with very deep X-ray observations. This combination allows 
us to use the full depth of the X-ray data, and to correct 
accurately for incompletenesses both in the X-ray detection 
and in the optical photometric selection. As a result, our 
determination of the faint end of the XLF should be both 
the deepest, and the most accurate yet obtained. Contrary 
to some previous results, and LDDE models, we find no 
evidence for a flattening of the faint end slope at high z. 
There is clear evolution of the XLF, however. This can be 
described in terms of zero or very mild positive luminosity 
evolution above L, , and negative density evolution below. 

At the very lowest luminosities probed by our survey 
(Lx < 10^^'^erg s~^) we have identifled a population of pure 
starburst galaxies, which are a significant contaminant of our 
XLF of AGN. Comparing with previous stacking results con- 
firms the existence of such a population, which is not found 
at z ~ 1, indicating a rapid increase in the number density 
of the brightest star-forming galaxies to high redshift. 

Our results are higher than previous estimates at the 
faint end, which we attribute to incompleteness in purely X- 
ray selected work. Our Lyman-break method identifies blue 
objects which, after correcting for incompleteness, accounts 
for all know X-ray sources at z — 3. This tentatively suggests 
that the AGN host population is evolving from blue star- 
forming galaxies at z = 3, to mainly red galaxies at z = 1. 

Future studies, utilising recent and forthcoming large 
area, deep surveys, will provide further constraints on the 
evolution of the XLF of AGN at high redshift. We stress 
that full consideration of selection and completeness effects 
is needed for a robust determination of the XLF. 
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